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A combined analysis of 3D Weak Lensing, Lyman-a forest and WMAP year three data 
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We present constraints on the amplitude and shape of the matter power spectrum and the den- 
sity of dark matter within the framework of a standard ACDM model. We use a Markov Chain 
Monte Carlo approach to combine independent measurements of the three dimensional weak grav- 
itational lensing shear field by the COSMOS survey, of low and high resolution Lya forest fiux 
power spectrum by SDSS and LUQAS, and of Cosmic Microwave Background temperature and 
polarization anisotropies by WMAP. We note good agreement between the amplitude of the matter 
power spectrum on intermediate and small scales as inferred from Lya forest and lensing data. The 
Lya forest data helps to break the erg — f^Om degeneracy characteristic of weak lensing results, yield- 
ing o-g = 0.876 ± 0.048 for COSMOS plus Lya SDSS data. This is somewhat larger than the value 
preferred by the WMAP year three CMB data. Combining all three data sets significantly tightens 
the constraints on erg, the spectral index of primordial density fluctuation rig, a possible running of 
the spectral index nrun and the matter density r^om. Assuming no running, the joint constraints for 
COSMOS, SDSS and WMAP are erg = 0.800 ± 0.023, = 0.971 ± 0.011, fiom = 0.247 ± 0.016 [l-a 
error bars). 

PACS numbers: 



I. INTRODUCTION 

The recent measurement of the Cosmic Microwave 
Background (CMB) anisotropies obtained by the WMAP 
satellite 0, [3 has considerably tightened the error 
bars on the cosmological parameters that describe the 
standard ACDM model of structure formation. The con- 
straining power of the WMAP year three measurements 
alone is already large. Cosmological parameter extrac- 
tion nevertheless benefits from a combination of the CMB 
data with measurements of the matter power spectrum 
on smaller scales and at different redshifts. The vari- 
ous observables suffer from very different systematic and 
statistical errors, and a combined analysis can break de- 
generacies inherent to individual measurements. 

We concentrate here on comparing and combining the 
CMB data with weak gravitational lensing and Lya forest 
data. The Lya forest due to the absorption produced 
by neutral hydrogen along the line-of-sight to distant 
quasars (QSOs) allows us to measure the matter power 
spectrum on scales that range from few to tens of co- 
moving Mpc at z = 2 — 4 (e.g. 0). "Cosmic shear" 
measurements of the distortion induced in distant galaxy 
images by weak gravitational lensing around foreground 
structures map the mass distribution at similar scales 
but smaller redshift, z = — 1.5. Both methods probe 
the matter power spectrum on smaller scales than CMB 
anisotropies observed by WMAP, and offer a more direct 
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measurement of the quantity erg, the r.m.s. of the den- 
sity fluctuations in spheres of radii 8 comoving Mpc. 
Furthermore, the small scale matter power spectrum in- 
ferred from the Lya forest data puts strong limits on the 
mass of warm dark matter candidates and on isocurva- 
ture perturbations (e.g. P, [13, [H, [10, HH ) . 

Vicl et al. [l^ and Seljak et al. [1^ have recently pre- 
sented independent analyses of CMB data combined with 
Lya forest data from separate, state-of-the-art QSO sam- 
ples: the LUQAS saniple of high-resolution high quality 
VLT-UVES spectra ( [ill, [H, [ll]) and a large sample of 
SDSS QSO spectra (McDonald et al. (H). Despite the 
very different data sets and the use of different analysis 
techniques, both groups found similar results, suggest- 
ing a value for cts ^ 0.9 larger than that extrapolated 
from the CMB data alone (erg = 0.76 ± 0.05). For the 
case of the SDSS Lya forest data there appears to be 
a moderate tension between the two data sets at the 2 
a level [10|. For the LUQAS data the errors are about 
a factor two larger and the difference is not statistically 
significant [l^. Ref. [lO] obtained very similar results 
with a further Lya forest data set. 

The moderate but notable tension between the am- 
plitude of the matter power spectrum inferred from the 
SDSS Lya forest data with the WMAP three year results 
suggests that measurements of the amplitude of the small 
and intermediate scale matter power spectrum are still 
somewhat uncertain. Further comparison with a third, 
completely independent technique should be very useful 
in this context. We exploit the tomographic weak grav- 
itational lensing analysis used to map the three dimen- 
sional distribution of mass in the Hubbble Space Telescope 
COSMOS survey [3i,[3^. The clustering signal in that 
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distribution also corresponds to a relatively large value of 
CTg = 0.951q (when the matter fraction flom is chosen 
to be the WMAP maximum likelihood value). A recent 
combination of two dimensional, ground-based cosmic 
shear surveys into the 100 square degree Weak Lensing 
Survey [3| gives a comparable although slightly smaller 
amplitude of the matter power spectrum erg = 0.84 ±0.07 
(for the best fit f^om WMAP value). 

In this paper, we investigate in detail to what extent 
the Lya forest and the cosmic shear data of Massey et al. 
(2007) are consistent, and present new constraints on the 
value of cosmological parameters obtained from various 
combinations of cosmic shear, Lya forest and CMB data 
(for the latter we limit ourselves to the WMAP year three 
determination of temperature/polarization anisotropics 
[III). We use the Boltzmann code CAMB [40] for comput- 
ing the linear matter power spectrum and applying rel- 
evant non linear corrections [39|. The multi dimensional 
parameter space is explored with Monte Carlo Markov 
Chains (MCMC) using the public code COSMOMC [H]. 



II. THE DATA SETS 
A. Lya forest data 

We have used two Lya forest data sets: i) the high 
resolution QSO absorption spectra presented in Viel, 
Haehnelt & Springel (VHS) ^ and in Refs. [H, [H, 
consisting of the LUQAS sample [ll[ and the reanalyzed 
Croft et al. (C02) data; ii) the SDSS Lya forest sam- 
ple presented in McDonald et al. [H (M05). The SDSS 
Lya forest data set consists of 3035 QSO spectra with 
low resolution {R ~ 2000) and low S/N (< 10 per pixel) 
spanning a wide range of redshifts {z = 2.2 — 4.2), while 
the LUQAS and the C02 samples contain mainly 57 high 
resolution {R ^ 45000), high signal-to-noise (> 50 per 
pixel) QSO spectra with median redshifts of z = 2.125 
and z = 2.72, respectively. 

The flux power spectrum of the Lya forest is the 
quantity which is observed and needs to be modeled at 
the percent or sub-percent level using accurate numeri- 
cal simulations that incorporate the relevant cosmologi- 
cal and astrophysical processes. M05 modeled the flux 
power spectrum using a large number of Hydro Parti- 
cle Mesh (HPM) simulations 0, ^ , calibrated with full 
hydrodynamical simulations. Instead, the VHS analysis 
significantly improved the effective bias method devel- 
oped by C02 (see Gnedin & Hamilton Q and Zaldar- 
riaga et al. [23] for a critical assessment of the errors 
involved) , by using a grid of full hydrodynamical simula- 
tions run with the Tree-SPH code GADGET-2 [H, [13] to 
infer the linear matter power spectrum. Finally, Viel & 
Haehnelt [2l| performed an independent analysis of the 
SDSS Lya forest data, and used a Taylor expansion of 
the flux power spectrum around best fitting values based 
on full hydrodynamical simulations to model the depen- 
dence of the flux power on cosmological and astrophys- 



ical parameters. This analysis was performed directly 
on the flux power spectrum and took thus full advan- 
tage of each data points, for a wider range of parameters 
than just the amplitude and slope constrained by the 
SDSS analysis. One should however keep in mind that 
the Taylor expansion approach is likely to underestimate 
the errors far from best-fit values. However, as soon as 
the SDSS Lyman-a data is combined with either COS- 
MOS or WMAP, large departures from the best-fit model 
are forbidden and the Taylor method is accurate. 

In this paper, we will use either the data set of 
VHS, based on high-resolution QSO spectra and noted 
"LyaVHS" ; or the results of Viel & Haehnelt [H us- 
ing low-resolution SDSS spectra and a Taylor expansion, 
noted "LyaSDSS-d" (where -d refers to "derivatives", 
since this method is based on the derivatives of the flux 
power spectrum. 

Our COSMOMC module lya. f 90 comparing the linear 
dark matter power with LyaVHS data has been incorpo- 
rated into the latest public available version of COSMOMC 
[28t . The LyaVHS power spectrum consists of estimates 
of the linear dark matter power spectrum at nine values 
in the wavenumber space k at z = 2.125 and nine values 
at z = 2.72, in the range 0.003 < k (s/km)< 0.03. The 
estimate of the uncertainty of the overall amplitude of 
the matter power spectrum is 29%. This estimate takes 
into account possible systematic and statistical errors 
(see the relevant tables of VHS for a detailed discussion). 
The code assigns a Gaussian prior to the correspond- 
ing nuisance parameter and marginalize over it. For the 
LyaSDSS-d analysis, we used the COSMOMC module de- 
scribed in [43 ]. which involves 21 nuisance parameters 
characterizing a wide range of astrophysical and noise- 
related systematic uncertainties. In the following results 
these parameters are always marginalized out. 



B. Weak lensing data 

The Hubbble Space Telescope COSMOS survey cov- 
ers a contiguous area of 1.64 square degrees on the sky. 
In this high resolution, space-based data, the shapes of 
234,370 distant galaxies were measured, with a median 
F814W0iB-band magnitude of 24.6 

A crucial addition to this survey has been the acquisi- 
tion of ground-based imaging in 15 extra bands. Photo- 
metric redshift estimation for each galaxy allows a fully 
3D exploitation of the signal, in which the power spec- 
trum can be independently measured at different red- 
shifts and physical scales. In a 2D analysis, these would 
have been projected together, resulting in a loss of dis- 
criminatory power on cosmological parameters of a factor 
of 3 — 5 [Sj, . BPZ photometric redshift estimation 
software achieved 68% confidence limits of 0.03(1 -I- z) 
on each galaxy to z ~ 1.4 and Ifsiaw = 24 [STf . 
The galaxy catalog has a median photometric redshift 
of ^phot = 1.26, and has been split into three redshift 
bins for this analysis: Zphot = 0.1 — 1, Zphot = 1 — 1.4 
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and Zphot = 1.4 — 3, which divide the number of sources 
almost evenly. 

In each redshift bin. the "cosmic shear" two-point cor- 
relation functions Ci,2(0) have been measured on angular 
scales 0.1—40 arcmin [33. Error estimates for these mea- 
surements include: statistical errors, calculated from the 
internal distribution of shear estimators; cosmic (sam- 
ple) variance, calculated from the variation in the sig- 
nal between separate quadrants of the COSMOS field; 
and finally systematic errors, which include the potential 
bias in shear calibration (overall bias and relative bias 
between redshift bins), errors due to catastrophic photo- 
metric redshift failures, and errors due to binning. 

In Massey et al. [sHi, the constraints on cosmolog- 
ical parameters were derived as follows. For a three- 
dimensional grid of models spanning variations of Qom 
from 0.05 to 1.1, erg from 0.35 to 1.4 and the power spec- 
trum shape parameter F from 0.13 to 0.33, the linear 
power spectra were obtained from the fitting formula of 
BBKS [4^, and corrected for non-linear evolution us- 
ing HALOFIT [3^. For each model, the data likehhood 
was computed taking only statistical errors into account. 
Then, the three-dimensional likelihood distribution was 
integrated in order to marginalize over F and to ob- 
tain confidence contours for (rJom, cs). The best con- 
strained combination of these parameters was found to 
be (78(riom/0.3)*' '*^. Final bounds on this quantity were 
obtained by adding systematic error linearly. 

In this work, we wrote a COSMOMC mod- 
ule for COSMOS data (downloadable at 
^ttp : //www. astro . caltech. edu/~rjm/ cosmos/cosmomc 
which offers various advantages with respect to the orig- 
inal analysis: the linear power spectrum is computed by 
CAMB, more free parameter can easily be implemented 
in the analysis (like the tilt ng), and the systematic 
errors can be accounted more accurately by introducing 
various nuisance parameters. The module computes the 
shear correlation functions for any cosmological model 
explored by the Markov chains (including corrections to 
the linear power spectrum obtained with the halofit 
code [13]) and compares with the tomographic results 
from Ref. [sHi, neglecting the cross-correlation between 
redshift slices. As shown in figure 8 of Ref. [sH], there 
are measurements on six angular scales, in each of three 
redshift bins. Since there are two correlation functions 
available for shear, the number of data points sums up 
to 36. The 36 x 36 covariance matrix is presented in 
figure 9 of [33 , and we used its inverse when calculating 
the likelihood. We modeled the systematics errors of 
the Icnsing data described in detail in Ref. [IHl by 
three nuisance parameters, over which our final results 
are marginalized. A blind analysis of simulated HST 
images suggests a potential 6% uncertainty in the overall 
calibration of the shear measurement |34l |. for which 
we account with a parameter A with Gaussian prior. 
We further introduce a parameter B to account for a 
5% relative calibration uncertainty between the shear 
measured from galaxies in our high and low redshift bins. 



Although not seen in the simulated HST simulations of 
[s^ l. more comprehensive tests on a larger set of simu- 
lated ground-based images [1, 01 reveal the potential for 
the shear to be underestimated in faint or small galaxies. 
The opposite effect has not been recorded. B is thus 
assigned a one-sided Gaussian prior. Finally, a potential 
10% intrusion of low redshift galaxies into the high 
redshift bin due to catastrophic photometric redshift 
errors is modeled with a third nuisance parameter C. 
Priors on the photometric redshifts were designed to 
ensure that the effect on lensing observables of known 
types of catastrophic failure are easily modeled. The 
failures can only dilute the signal in the high redshift 
bin, and again we use a one-sided Gaussian prior. In 
summary, the data points Ci.2(0) are multiplied by 

{A/By 
A^ 



[ABf/C 



(low— z) 
(medium- z) 
(high-z) 



with priors on A, B and C peaking at one with S > 1, 
C>l,cjA = 0.06, (IB = 0.05 and ac = 0.10. 



III. RESULTS 

For our MCMC analysis we have assumed a minimal 
flat ACDM model, with no tensor contribution. We 
vary the following cosmological parameters with top- 
hat priors: dark matter density ^och^ G [0.01,0.99], 
baryon density flobh^ G [0.005,0.1], primordial spectral 
/)index rig 6 [0.5, 1.5], primordial amplitude log[10^°As] G 
[2.7,4.0] and angular diameter of the sound horizon at 
last scattering 9 G [0.5,10]. When CMB data is in- 
cluded, we also vary the optical depth to rcionization 
r e [0.01,0.8]. For part of the MCMC analysis, we 
have finally considered a running of the spectral in- 
dex, riiun e [—0.5,0.5]. From the parameters above 
the MCMC code derives the reduced Hubble parameter 
h, the matter fraction f2om and erg: so, these parame- 
ters have non-flat priors and the corresponding bounds 
must be interpreted with some care. In addition, COS- 
MOMC imposes a weak prior on the Hubble parameter: 
h G [0.4, 1.0]. In each case and for the purpose of com- 
paring with Ref. (35| . we also compute the combination 
CT8(f^om/0.3)°-4'', which is best probed by the COSMOS 
lensing data. 

We ran first a MCMC analysis of the COSMOS WL 
data alone. We recall that this analysis differs from 
that in Massey et al. [IH] in four ways: the method- 
ology (MCMC with flat priors on the cosmological pa- 
rameters instead of maximum likelihood), the introduc- 
tion of one extra independent parameter rig, the fact 
that we compute the exact linear matter power spec- 
tra numericall y, and the treatment of systematic er- 
rors. In Ref. [35|, the analysis performed with statis- 
tical errors only led to the bounds erg (^^Om/0. 3)" "''^ — 
0.866 ±0.033 (68% Confidence Level (C.L.)), which ex- 
tend to cr8(f^0m/0.3)" '*'^ = 0.866 j:°;o68 when systematic 
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Figure 1: 1-a and 2-a contours of the marginalized likelihood in the ag — flom plan e (left) and erg — plane (right) for each 
data set separately: WMAP (red), weak lensing (orange), Lya forest from VHS [22{ (blue) and from SDSS as analyzed by [2ll | 
(green). 




Figure 2; 1-a and 2-a contours of the marginalized likelihood in the erg — Jlom plane (left column) and as — ria pi 
column) for various combinations of data sets from COSMOS (WL), WMAP, and Lya forest data from VHS [22| 
upper panels) and SDSS as analyzed by [2l| (LyaSDSS-d, lower panels). 



ilane (right 
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Table I: Summary of the constraints on erg, Us, Slom, h and 
r, for the minimal 6-parameter ACDM model and each data 
sets. Since this is a Bayesian analysis, the bounds depend on 
our choice of priors; our top-hat priors are described at the 
beginning of Sec. lIIII fin particular, we impose a weak h prior: 
0.4 < h < 1.0). The quoted values are either the mean and 
68% C.L. error, or only the 68% C.L. upper/lower limit when 
a parameter is not bounded on both sides within the prior 
range. 





WL 


Lyn VHS 


Lya SDSS-d 


WMAP3 




0.85 ±0.22 


1.04 ±0.16 


0.926 ±0.066 


0.762 ±0.046 


Tls 


< 0.94 


0.80 ±0.10 


0.982 ± 0.028 


0.955 ±0.016 




0.34 ±0.19 


0.55 ±0.26 


0.238 ± 0.030 


0.243 ± 0.032 


h 


> 0.71 


> 0.63 


0.710 ±0.071 


0.729 ±0.030 


T 








< 0.104 



Table II: Same as table |I] for the combination of weak lensing 
data with each other data set. 



WL+WMAP3 WL + Lya VHS WL + Lya SDSS-d 



OS 


0.802 ± 0.034 


0.98 ±0.19 


0.876 ± 0.048 


Ha 


0.958 ±0.016 


0.88 ±0.11 


0.962 ± 0.034 




0.269 ±0.026 


0.25 ±0.12 


0.232 ± 0.028 


h 


0.708 ± 0.023 


> 0.79 


0.773 ± 0.089 


r 


< 0.103 







errors are added linearly. Here wc first performed a 
MCMC analysis excluding the systematic errors, and ob- 
tained cr8(^^0m/0. 3)°-'"' = 0.83±0.06. The increased error 
bar with respect to [sst is presumably due to the inclusion 
of an arbitrary spectral index rig in the analysis, which 
opens new parameter degeneracies. The mean value of 
Massey et al. (2007) (0.866) is perfectly consistent with 
this value. We then incorporated systematic errors as 
described in the last section and found 

erg (f^om/O.S)" "*^ = 0.814 ±0.074 (68%C.L.) . (IILl) 

Note that because of non-linear corrections to the matter 
power spectrum, the erg parameter cannot be viewed as 
a simple calibration parameter for the theoretical corre- 
lation functions. A change in trg changes both the ampli- 
tude and the shape of the shear correlation functions in a 
non-trivial way. As a consequence, the impact of system- 
atic errors on the determination of the (linear theory) 
parameter (T8(r^om/0.3)" '*"' is found to be smaller than 
the data calibration uncertainty itself. The 68% confi- 
dence limits on each parameter arc presented in the first 
column of Tabic I. Wc also show the joint 68% and 95% 
confidence contours in crg-fiom and ag-ris space in Fig. [1] 
(yellow). 

For comparison, we also show in Fig.[T]the as-i^om and 
ag-Hg confidence regions for each dataset separately. The 
COSMOS data is compatible with WMAP, since the con- 
tours have some overlap even at the 68% level. In the a^- 
fiom space, the COSMOS and WMAP contours appear as 



Table III: Same as table|T]for the combination of CMB, weak 
lensing and Lya forest data. The quoted values are the mean 
and 68% confidence limits. 



WL + WMAP3 + Lya VHS WL + WMAP.S + Lya SDSS-d 



0"8 


0.822 ±0.032 


0.800 ± 0.023 




0.960 ± 0.016 


0.971 ±0.011 




0.282 ± 0.026 


0.247 ±0.016 


h 


0.700 ± 0.022 


0.730 ± 0.016 


r 


0.094 ± 0.028 


0.109 ±0.026 



Table IV: Same as table Hill for the case with a running spec- 
tral index. The quoted values are the mean and 68% confi- 
dence limits. 



WL + WMAP3 + LyQ VHS WL+WMAP3 + Lyn SDSS-d 



0"8 


0.809 ± 0.041 


0.818 ±0.024 


ris 


0.965 ± 0.018 


0.971 ±0.015 




0.304 ± 0.032 


0.255 ± 0.018 


h 


0.679 ± 0.026 


0.719 ±0.018 


r 


0.085 ± 0.037 


0.135 ±0.026 




-0.028 ±0.018 


-0.007 ±0.021 



almost orthogonal; and the overlap clearly suggests that 
the WL data prefers the highest erg values allowed by 
WMAP. The COSMOS data is also compatible with the 
different Lya data sets, with again an overlap at the 68% 
level. Finally, the WMAP and Lya contours overlap only 
at the 2 a level, as expected from previous works [2^, \^ . 

The two panels in Fig. [T] provide a good illustration of 
the advantages of combining various datasets. For each 
type of experiment, erg and Oom are clearly correlated, 
but the direction of the degeneracy is different for CMB, 
WL and Lya data. There are various reasons for the 
difference between the direction of correlation associated 
with the WL and Lya data: first, the raw Lya data are 
in units of s/km (since the power spectrum is measured 
in velocity space), and the rescaling to units of h/Mpc 
depends on fiom; second, the data probe different red- 
shifts, and the ratio between the power spectrum today 
(when as is defined) and at a given redshift depends on 
f2om; third, the slope of the matter power spectrum de- 
pends on fiom • This explains the "banana shapes" in the 
upper left panels, with different orientations. 

We don't find any significant correlation between erg 
and Us (upper right panel in Fig. [1]). For the WL and 
Lya data, this is due to the fact that these experiments 
directly probe power on the scale at which cg is defined 
(if this was not the case, the amplitude of the WL and 
Lya experimental points would constrain a combination 
of (Tg and Us). 

Wc then ran some Markov chains for various combi- 
nations of different data sets and we show the results in 
Table [U IIIII and Fig. [2] Combining the data sets sig- 
nificantly tightens the constraints. Most noteworthy is 
that there remains a (moderate) strain between the in- 
ferred value of (Tg between that inferred from WMAP 
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alone and that inferred from the lensing and Lya forest 
data. The constraint on ag from COSMOS+LyaVHS 
(erg = 0.98 ± 0.19) and from COSMOS+LyaSDSS-d 
{as = 0.876 ± 0.048) are compatible with the WMAP 
best-fit value (erg = 0.762) respectively at the l.l-cr and 
2.4-(T level. Note that the SDSS-d contours are based 
on extrapolations using a Taylor expansion of the flux 
power spectrum around a best fit model and arc likely 
to underestimate the error for parameters far from the 
best fit model. However, when this data set is used 
in combination with WL and/or WMAP data, the 68% 
and 95% C.L. contours remain in a small region where 
the Taylor expansion is accurate. The main results of 
this work are the 68% confidence limits for the com- 
bined analysis of CMB, weak lensing and Lya forest 
data: WMAP-hCOSMOS-hLyaVHS, erg = 0.822 ±0.032; 
WMAP-FCOSMOS-hLyaSDSS-d, ag = 0.800 ± 0.023. 

Finally, we performed a further MCMC analysis with 
an extended ACDM model with on extra parameter, a 
running of the spectral index nj-un- The results are sum- 
marized in Table IV. In this case, the tilt Ug is defined at 
the pivot scale fco = 0.01 Mpc~^ (when WL and Lya data 
are included, the pivot value fco = 0.002 Mpc~^ adopted 
in the WMAP3 paper [31 is too small with respect to 
the median scale of the full data set). The choice of a 
given pivot scale is indifferent for the definition of Uj-un- 
WMAP alone is compatible with a rather large negative 
running, rijun ~ —0.055 ± 0.03 (68%C.L.) which results 
in a reduction of power on small scale. Small scale exper- 
iments like Lya forest and weak lensing observations are 
obviously crucial for the determination of a possible run- 
ning of the spectral index, since they increase the lever 
arm for primordial spectrum reconstruction. The high erg 
value preferred by WL and Lya data excludes the most 
negative values of rirun found by WMAP, and reduce its 
error by a factor of two. Wefindrii un = — 0.028±0.018 for 
WMAP+COSMOS+LyaVHS and n^n = -0.007±0.021 
for WMAP+COSMOS-l-LyaSDSS-d (68%C.L.). 



flux derivatives method of Vicl & Hachnelt [2l|) than 
the WMAP data alone (erg = 0.762 ± 0.046). The 
direction of degeneracy between the amplitude of the 
power spectrum on galaxy scales the parameters gov- 
erning its shape (in other words, the direction of degen- 
eracy in cTg-ng and crg-f2om space) is different for the 
Lya forest and weak lensing data. These two obscrv- 
ables thus complement each other very well and combin- 
ing them results in a significant improvement. Combin- 
ing all three observables we get either as = 0.800 ±0.023, 
71s = 0.971 ±0.011, = 0.247 ±0.016 (with the analy- 
sis of SDSS Lya forest data based on the flux derivatives 
method of Viel & Haehnelt [2l|) or as = 0.822 ± 0.032, 
ris = 0.960 ± 0.016, Qom = 0.282 ± 0.026 (with the high 
resolution Lya data of VHS). We further explored the 
constraints for a running of the spectral index and found 
the data to be consistent with no running at less than 
the 2a level. Adding the smaller scale data sets reduces 
the uncertainty on the running of the spectral index by 
a factor of two with respect to WMAP alone. 
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IV. CONCLUSIONS 



We have presented a joint analysis of the constraints on 
the matter power spectrum and the density of dark mat- 
ter from three different cosmological probes: the CMB 
temperature and polarization anisotropy measurements 
of WMAP year three, the state-of-the-art weak lensing 
COSMOS survey and two independent Lya forest data 
sets. The different observables are prone to very different 
systematic errors and parameter degeneracies, and more 
importantly probe different scales and rcdshifts. Assess- 
ing their consistency is an important test of the ACDM 
paradigm and is crucial for further improving constraints 
on cosmological parameters. The measurements of the 
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